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ABSTRACT 



Context. Recent attention has been directed to abundance variations among very young stars. 

Aims. To perform a detailed abundance study of the Herbig Ae star HD 101412, taking advantage of its unusually sharp 
spectral lines. 

Methods. High-resolution spectra are measured for accurate wavelengths and equivalent widths. Balmer-line fits and 
' ionization equlibria give a relation between Teg, and log(g). Abundance anomalies and uncertain reddening preclude 

C-j ' the use of spectral type or photometry to fix T e g. Excitation temperatures are used to break the degeneracy between 

O ,1 T off and log(g). 

Results. Strong lines are subject to an anomalous saturation that cannot be removed by assuming a low microturbulence. 
By restricting the analysis to weak (< 20mA) lines, we find consistent results for neutral and ionized species, based 
on a model with T e g = 8300K, and log(g) = 3.8. The photosphere is depleted in the most refractory elements, while 
volatiles are normal or, in the case of nitrogen, overabundant with respect to the sun. The anomalies are unlike those 
of Ap or Am stars. 

Conclusions. We suggest the anomalous saturation of strong lines arises from heating of the upper atmospheric layers 
by infalling material from a disk. The overall abundance pattern may be related to those found for the A Boo stars, 
though the depletions of the refractory elements are milder, more like those of Vega. However, the intermediate volatile 
zinc is depleted, precluding a straightforward interpretation of the abundance pattern in terms of gas-grain separation. 

Key words. Stars: chemically peculiar - Stars: abundances - Stars: individual (HD101412)- Stars:pre-main-sequence 

0© 1. Introduction ally sharp, leading to suggestions that the star may be seen 

< 5 mn lm ,^,t^ ™ nn r,r i r- n m ^ \ , , nearly "pole on." Intrinsically slow rotation cannot be ex- 

. The star HD 101412 ! (CD -59 3865, V1052 Cen) belongs ^ ^ Fedd gt &l y {2QQ8) & ^ ^ 

r-H . to the group of Herbig Ae/Be stars, whose members are & digk b 2QO to the Une of gi M a model now 

> considered more massive counterparts of T Tauri pre-mam supported b a new l y -fo U nd photometric period. 

sequence stars. Masses range from 2-10 Alp, (Hernandez, ct ,, , , , . JTr . ,„,,,„ . , 

1 j, , n ™r\ mi • i l- i -i y ■ i 1 Very recent work shows that HD 101412 is a low- 
^ ■ al. 2005). Iheir spectral energy distribution is character- ,., J , . . . , , ... . , c Anrir , 
r ^i r . , ■ , j .i i amplitude photometric variable with the period ol 42?0 

I— ( lzed by the presence of an infrared excess due to thermal , , . , , , . , . N Tr , , . , . 

rrt ■ . . c ■ ,11 i i • i • ,i , , , (Mikulasek et al., in preparation). If the period is a ro- 

w ^ , re-emission ol circumstellar dust, which is thought to be > , , ., -,i 

, c , „ j- i / r TT i onnn tational one, it then comports well with our spectroscopic 

the signature of a circumstellar disk (cf. Hartmann 2009, , . . „,.,•, _i , , , , . 

,oo\a j- nji ,i )nnna\ Tm mi ^10 • value, v • sin i = 3 ± 1 km s , and excludes a conclusion 

beet. 8.8 . According to Fedele et al. (2008 HD 101412 is , , , ' ,. , , , ' , , / , , , 

' ,. , .„ \. r at i that the narrow lines could only be due to a (nearly) pole-on 

a group II source (using the classification of Meeus et al. , v ' 

(2001), with the flat disk self-shadowed by dust. Wade, et V1CT ™§ an S le - . . „ . TT , . 4 

al. (2005) discussed HD 101412 as a possible progenitor of Numerous emission features m Herbig Ae stars indicate 

magnetic chemically peculiar (CP2 or Ap) stars. A number an , atmospheric structure that must differ from the classi- 

of recent studies have been devoted to this star (cf . Hubrig, cal ' , one-dimensional models typically used for abundance 

et al. 2009, 2010), including an estimate of its metaUicity f nal y ses - Nevertheless abundance work on these objects 

(Guimaraes, et al 2006). The absorption lines are unusu- 11 been a A t ^ n P t ^ d (<£ Guimaraes, et al 2006, Acke and 

Welkens 2004). I he latter study is of particular interest to 

Send offprint requests to: C. R. Cowley the results presented here, as the authors sought the sig- 

* Based on observations obtained at the European Southern nature of the A Boo phenomena, citing the seminal paper 

Observatory, Paranal and La Silla, Chile (ESO programmes by Venn and Lambert (1990). These authors noted that the 

077.C-0521(A), 081.C-0410(A) and 383.C-0684(A) most volatile elements, C, N, and O showed mild abundance 
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depletions, while the more refractory elements showed sig- 
nificant depletions, 1 dex or more. They noted the simi- 
larity of these depletions to those of the interstellar gas 
which is well correlated with condensation temperatures 
(cf. Sect. [7]). 

The sharp spectral lines of HD 101412 make it ideal 
for abundance work. Although a magnetic field has been 
detected (cf. Wade, et al. 2005), magnetic enhancement of 
the lines is masked by other factors, to be discussed below 
(Sect.[S]). Hubrig et al. (2010) noted that the Zeeman pat- 
terns of several lines were resolved, or partially resolved. 
Nevertheless, most atomic lines are sufficiently sharp that 
numerous relatively unblended lines could be found for 
abundance work. 

2. Spectra 

The current paper is based on a subset of the spectra de- 
scribed by Hubrig et al. (2010). Most equivalent widths 
were measured on two HARPS spectra (averaged) obtained 
in 2009 on 23:47 UT of 4 July, and 00:13 UT of 5 July. 
About 10 per cent of the equivalent widths are from UVES 
spectra obtained on 9 (visual and IR) and 15 April (UV) 
2009. 

Wavelengths were measured on the HARPS spectrum 
in the range AA3782-6911. The averaged spectra were sub- 
jected to mild Fourier filtering. Signal-to-noise estimates, 
made directly from the filtered spectra average 148. HARPS 
spectra have a resolving power (RP) of 120000, but av- 
eraging and filtering reduce this to an estimated 40000. 
Fortunately, equivalent widths are independent of RP; we 
exchange RP for noise reduction. 

We have used accurate wavelength measurements both 
for line identification, and to find relatively unblended lines. 

The UVES spectrum obtained on 15 April 2009, was 
measured for wavelengths in the range AA3301-4517. A few 
equivalent widths below the Balmer Limit were obtained 
from this spectrum. 

Additional wavelengths were measured from UVES 
spectra taken on 9 April 1999 (AA6911-7519 and 7661- 
9461). The Ca n infrared triplet was examined for possible 
isotopic shifts discovered by Castelli and Hubrig (2004). 
The lines were shifted by 0.02 to 0.03 A with respect to 
the solar wavelengths of these lines, a rather small shift 
that could be entirely due to measurement error. 

Hubrig, et. al. (2010) demonstrate equivalent width 
changes of the order of 20% (0.08 dex). The variations were 
correlated with an assumed phase, 13^86 days nearly one 
third the newly determined period of 42*?0. Also possible, 
are line strength changes due to time-dependent motions of 
the circumstellar material. These variations must be kept in 
mind in evaluating the abundances which were determined 
here primarily from a single phase. An error of 0.08 dex is 
marginally significant, relative to other uncertainties in the 
abundances (cf. Tabled]). 

Balmer line profiles were measured on a FORS 1 spec- 
trum obtained on 22 May 2008, in visitor mode at ESO. We 
used the (200kHz, low, lxl) readout mode, which makes 
it possible to achieve a S/N ratio of about 1000-1200 with 
a single exposure and the GRISM600B in the wavelength 
range 3250-6215 A to cover all hydrogen Balmer lines from 
H/3 to the Balmer jump. A slit width of 0'.'4 was used to 
obtain a spectral resolving power of R s=s 2000. Details of 
data reduction are given by Hubrig et al. (2004). 



3. Model atmospheres and spectral calculations 

All model atmospheres are based on T—t relations from the 
version of ATLAS9 (Kurucz 1993) of the Trieste group (cf. 
Sbordone, et al. 2004). Given T(r), the depth integrations 
from log(r5ooo) of —5.4 to +1.4 were carried out with soft- 
ware and opacity routines written and used at Michigan for 
several decades. Agreement with ATLAS9 Models posted 
on Castelli's (2010) site are excellent (see Castelli & Kurucz 
2003). Spot checks using WIDTH9 for weak Fe I and n lines, 
which are independent of the different damping constants 
of the two codes, show agreement within 0.02 dex. Weak 
lines have equivalent widths of 20 mA or less. 



4. Fixing the atmospheric parameters 

4.1. Colors and spectral type 

Broad and intermediate-band photometric measurements 
are typically used to fix the effective temperatures prior to 
an abundance analysis. Johnson (Vieira, et al. 2003) and 
Geneva (Mermillod, Hauck, and Mermillod 2007) indices 
are available for HD 101412; we found no Stromgren pho- 
tometry. Unfortunately, the reddening correction is most 
uncertain. For example, Vieira, et al.'s U — B = 0.15, 
B — V = 0.18 do not lead to an intersection with a standard 
U — B vs. B — V plot for either main sequence stars (Cox 
1995, Table 15.7) or a ZAMS (Table 15.9). Since the star is 
known to be peculiar and variable, currently available pho- 
tometry is not adequate to fix the effective temperature. 

The spectral type of HD 101412 is given variously as B9- 
A0 V-III (e.g. B9/A0 V, Houk and Cowley 1975). According 
to Cox (1995), luminosity class V stars have effective tem- 
peratures in the range 10500 to 9790K, while we conclude 
T e ff ?» 8300K. The spectral type of normal stars of this class 
depend strongly on the Ca II K-line, or rather its strength 
relative to Ca n H +He. The intensity of the Ca n lines 
would be significantly diminished by the « 0.6 dex under- 
abundance we find for this element. 

We conclude that the spectral type is not a useful indi- 
cator of the temperature of the star, or of the total absorp- 
tion, Ay. 

4.2. The Balmer lines 

Balmer line profiles are difficult to work with on high- 
dispersion material. This is especially true for echelle spec- 
tra, as the broad lines typically span echelle orders. We 
have attempted to fit profiles from an averaged HARPS 
spectrum, but have relied primarily on the FORS1 spec- 
trum (Sect. [J) which included H/3 and higher members of 
the Balmer series (Fig. [1]). 

Guimaraes, et al. (2006) made use of Balmer lines but 
did not discuss T e g — log(g) degeneracy. They conclude 
T cS = 10000 ± 1000K (see their Table 2). In Hubrig, et 
al. (2009), only the wings of H/3 were used to find the at- 
mospheric parameters, T e g = 10000, log (g) — 4.2—4.3. The 
authors note that a 9000K, log (g) — 4.0 model also gave a 
good fit to the H/3 profile, but that with these parameters 
"many narrow lines" were seen in the synthetic spectrum 
that did not appear in the observed spectrum. These lines 
might not have appeared in the calculated spectrum if the 
abundances were low, as we currently assume. 
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Table 1. Provisional abundances (\og(Fe/N to t)) from Fe 
I and Fe II for various temperatures and surface gravities. 
Nearly equal abundances for Fe I and n are in bold face. 
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Fig. 2. Values of T c g/100 and log (g) yielding equal abun- 
dances for Fe I and Fe II. 



Wavelength [A] 

Fig. 1. H/3 (above, top wavelength scale) and Hi5 profiles 
(bottom wavelength scale) for T e g = 8300K, log(g) = 3.8 
and 9800K, log(<?) = 3.8. Observed profiles are in gray (red 
in online version) with dots. The fits for the lower temper- 
ature are as good or arguably better than at the higher 
temperature. The relative flux scale is correct for the low- 
est profile pair. Upper pairs are successively raised by 0.4 
units. 

The Hubrig et al. (2009) results show that within the 
relevant temperature range for HD 101412, the Balmer pro- 
files are sensitive to both temperature and gravity. Our cur- 
rent calculations show that quite good fits to the Balmer 
profiles may be obtained for a temperature as low as 8300K, 
if a log (g) = 3.8 is assumed (see Fig. [T]). A still lower 
temperature could be accommodated by assuming an even 
lower gravity. 

Wavelength measurements of the cores of the Balmer 
lines on the FORS1 spectra show systematic red shifts that 
increase systematically from about +9 at H19 to +23 km 
s _1 at H/3. Core shifts are marginally detectable in Fig. Q] 

4.3. Ionization equilibrium 

We assume here that an optimum model will yield the same 
abundances from two stages of ionization of a given ele- 
ment. In a preliminary approach to obtain model parame- 
ters, we adopted a grid of temperatures and gravities shown 



in Table [TJ The grid was chosen to bracket the expected 
stellar parameters. Provisional abundances were then cal- 
culated from equivalent width measurements of Fe I and n 
lines using models with the parameters given in the table. 

A relation between T e s and log (g) is obtained as follows. 
For each of the three values of log (g), we plot \og(Fe/N tot ) 
vs. T c ff and find a value of T e g where the abundances from 
Fe I and Fe II agree. This is illustrated in Fig. [2j 

5. Anomalous saturation 

The abundance from a set of equivalent widths of a given 
spectrum, e.g. Fe I, or Ti II, should not, in principle, de- 
pend on the line strength. However, it has been known since 
the early days of curve of growth analysis, that drifts of 
abundance with line strength are common. They can result 
from a variety of causes ranging from errors in the oscillator 
strengths to equivalent width measurements. In a common 
situation, strong lines will give a higher abundance than 
weak ones, and the abundance worker can remove the trend 
by assuming an additional source of broadening known as 
microturbulence (£t). This broadening arises from gas mo- 
tions due to convection, mixing, or other sources not in- 
cluded in the basic expression for the absorption profile, 
such as hyperfine structure or Zeeman effect. 

Plots of abundance vs. equivalent width for HD 101412 
slope downward even when no microturbulence is assumed. 
Fig. [3] shows the behavior for Ti II lines. Not all plots show 
the flattening for weak lines. In particular, the Fe II plot 
does not. However, most do. 
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Fig. 3. \og(Ti/N tot ) vs. log [W\{mA)} for Ti n lines using 
the model with T c g = 8300K, log(g) = 3.8. A microturbu- 
lence £t = 0, was assumed. The downward trend of abun- 
dance with equivalent width illustrates the anomalous sat- 
uration. (Additional plots may be found at the URL given 
in Sect. 

An obvious mechanism that would weaken lines in 
young stars is veiling. This arises from excess continuum 
radiation arising from some aspect star-formation process 
(cf. Hartmann 2009). It is important to assess the possible 
effects of veiling, which would weaken the absorption lines 
and mimic lower abundances. 

We believe the effect of veiling to be minor, for the fol- 
lowing reasons: 

• Carbon, oxygen, and sulfur give abundances close to 
solar. In the case of nitrogen, we find a significant ex- 
cess. Veiling, if it were present, should have affected lines 
from these elements. 

• Veiling is often wavelength dependent. There is no indi- 
cation of a drift in our abundances with wavelength for 
the five spectra having lines on either side of the Balmer 
jump: Mg I, V II, Cr I, Co I, and Co II. 

• Veiling is caused by an excess continuum due to infall 
from a disk. The veiling is less important for the Herbig 
Ae stars, relative to the T Tauri or weak T Tauri stars. 
There is no indication of X-rays from HD 101412, which 
would accompany significant infall of material. 

We have found that anomalous saturation does not oc- 
cur in a model with a hot upper atmosphere. If we arbitrar- 
ily raise the temperature from the standard ATLAS9 T(t) 
the effect vanishes (see Fig. 0]). Temperature changes were 
made by trial and error until the anomalous saturation was 
no longer evident in a plot of abundance from Fe II lines vs. 
equivalent width. The altered model is surely not unique. 
Since the results are from an LTE calculations, we take it 
as indicative of the problem with the ATLAS9 model for 
HD 101412 specifically. A realistic pursuit of this problem 
is best left to a NLTE study. 



6. T e g — log (<?) degeneracy; curves of growth 

We can see from Fig.[U and from Table [TJ that higher tem- 
peratures require higher gravities. This is true for both the 




lo S( T 500o) 

Fig. 4. Temperature distributions for the ATLAS9 model 
with T eS = 8300K, log (g) = 3.8 (solid line), and a "hot 
top" model (dashed) which gives no anomalous saturation 
for & = 0. 



iron equilibrium and the Balmer profiles. The temperature 
(9500-10000K) obtained by earlier workers would not yield 
acceptable Fe ionization, though the Balmer lines could be 
accommodated with a higher surface gravity. 

This leaves the question of how to break the degeneracy 
between temperature and gravity. In one of the relatively 
small number of abundance analyses of pre-main sequence 
stars, Acke and Waelkens (2004) determined atomic exci- 
tation temperatures "for ions with many observed lines." 
Specifically, they required that the effective temperature of 
the model be such that there is no drift in the abundance 
determined from individual lines with excitation potential. 

We have attempted to use the Acke and Waelkens 
(2004) method. Our results marginally favor a 8300K- 
log(g) = 3.8 model over a 9800K-log(#) = 4.2 model, 
though not decisively. The wrong temperature should show 
a clear difference between abundances from low- and high- 
excitation lines. Unfortunately anomalous saturation pre- 
cluded use of many stronger lines, which are typically of 
low excitation. As we shall see, this was not a problem with 
the second method, based on a curve of growth technique. 

Since our depth-dependent models are not accurate for 
this Herbig Ae star, it is worthwhile to explore results of 
a more basic approach where the photosphere is approx- 
imated by a uniform slab, with a single, mean value of 
the temperature and pressure. This simple technique was 
once widely used in the chemical analysis of stellar spectra, 
employing what were called Schuster-Schwarzschild (hence- 
forth, SS) models, (cf. Aller 1963). The basic method is still 
routinely employed when the conditions along the line of 
sight are not well determined (cf. Spitzer 1978, Rachford, 
et al. 2001, Hobbs 2005). This is the case for the upper 
layers of HD 101412. 

Note: we use the SS models only to find an excitation 
temperature and break the T e g — log(g) degeneracy. 

To determine an excitation temperature, we plot 
log (W A /A)+6 vs. log (gfX)-e-x- Here, 9 = 5040./T OXC , and 
X is the lower excitation potential in eV. One adjusts the 
value of 9 until the points belonging to lines with different 
excitation delineate, as well as possible, the same curve. 
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Constants are added to both the ordinates and abscissae 
of the observed points, in order to make the theoretical 
and stellar curves of growth overlap. The ordinate of the 
theoretical curve is \og\W\/ {2r A\n)\, while that of the 
empirical plot is log(WA/A) + 6. The vertical shift necessary 
to superimpose the two curves contains information on the 
Doppler width, A Ad = |\/ 27 ^ T /V + and the maximum 
line depth ro, discussed in more detail below. The horizontal 
shift contains information on the column density, but that 
is not needed here. The analytical curves were originally 
due to van der Held (1931), and are are shown in the figure 
as solid lines for several values of the ratio of a — 7a /(2 • 
AAd). Here 7a is the damping constant (in cm), /i is the 
molecular (atomic) weight, and 1Z the gas constant. Note 
that a constant has been added to our abscissa so that X 
and the ordinates of the theoretical curves are the same for 
weak lines. We made curves of growth for Ti II, Fe I, and Fe 
II, spectra with better-quality oscillator strengths, and Ca 
I We have tried to select oscillator strengths of optimum 
quality. Sources are: 

• Ti ii. Oscillator strengths are from Pickering, Thorne & 
Perez (2001) 

• Cr I. Oscillator strengths from Sobeck, Lawler, & 
Sneden (2007). 

• Cr II Oscillator strengths from Nilsson, et al. (2006) or 
VALD (Kupka, et al. 1999)/Kurucz (1994) taking only 
LS-permitted lines. Results were essentially the same. 

• Fe I: Lines from the NIST site (Ralchenko, et al. 
2010) with accuracy B+ (7 lines) and C+ (11 lines). 
Unfortunately, there are no lines with excitation poten- 
tials above 2.6 eV of accuracy B+ (see the NIST site 
for an explanation of "accuracy.") 

• Fe II: We adopt the oscillator strengths of Melendez and 
Barbuy (2009). 

Fig. [5] shows a curve of growth for Fe II. Too large a 9 
(too low a temperature) will move the filled stars (highest 
x) to the left with respect to the open circles (lowest x)- 
Too small a 9, would move the highest-x points to the right. 
Points for the lowest-excitation lines move relatively little. 

In the case of Ca I, there is no good recent set 
of oscillator strengths, apart from the three lines given 
by Aldenius, Lundberg, and Blackwell- Whitehead (2009), 
AA6162.2, 6122.2, and 6102.7). These three lines are all of 
moderate strength (9.6 to 25.3 mA), and give —6.41, —6.42, 
and —6.30 for \og(Ca/N to t)- We adopt —6.38; the mean of 
values from 35 lines is —6.43 ± 0.23 sd. A curve of growth 
is shown for Ca I in Fig. [5] Oscillator strengths for the re- 
maining 32 lines are from the NIST site (Ralchenko, et al. 
2010). 

Additional curves of growth and examples showing re- 
sults of a wrong temperature may be found at the url given 
in Sect. [3 

Curves of growth provide a measure of the excita- 
tion temperature, the microturbulence (£t), and the damp- 
ing constant. The excitation temperatures found from the 
curves of growth, corresponding to 9 = 0.65 to 0.7 (T = 
7750 to 7200K) are much cooler than any temperature dis- 
cussed so far for HD 101412. However, the excitation tem- 
perature from the line spectrum is expected to be signifi- 
cantly lower than that from the continuum. We have made 
calculations based on model atmospheres to verify that the 
above excitation temperatures are to be expected for an 
atmosphere with T e g = 8300K. 




X 

Fig. 5. Curve of growth for Fe II lines using oscillator 
strengths from Melendez and Barbuy (2009). Unit slope 
is indicated by the solid gray line. Along this line, the ordi- 
nate, and the abscissa {X), are equal. SS model curves are 
drawn for a = 7A /(2 • AA D ) of 1, 0.1, 0.01, and 0.001 (from 
the top). The abscissae for the points are log (g/A) — 9 ■ X- 
Open circles are for 2.58 < x < 2.89, open squares for 
3.20 < x < 3-89, and filled stars for 5.51 < x < 6.22 eV. 
The observed points have been made to coincide by adding 
+5.7 to the abscissa, and —0.6 to the ordinate. The lat- 
ter figure is relevant for the microturbulence. This plot was 
made with 9 = 0.7 (7200K). 




X 

Fig. 6. Curve of growth for 35 Ca I lines. The resonance 
line, A4227 is the open circle. Open squares are for the 
range 1.88 <x< 1-90, triangles 2.50 <x< 2 -71, and filled 
stars (2 lines) 2.93 eV. Ax = 1.7 Ay = -0.65 and 9 = 0.68 
(T = 7410K). Theoretical curves as in Fig. [5] 

The vertical displacement (Ay) necessary to superim- 
pose the theoretical and empirical curves gives a relation 
between the maximum line depth ro, an assumed relevant 
temperature, and the microturbulence. 

Our problem with the SS-model curves of growth is es- 
sentially the same as that with the ATLAS9 model. The 
maximum line depth expected for either model is too large, 
relative to the value needed in HD 101412 to give a finite 
microturbulence. In the case of the numerical model, ro in 
LTE is simply [1- B x (T min )/(2F§)}. The value, at A4481 in 
the 8300K- log(c?) = 3.8 model is 0.90. For the SS model, 
if we adopt ro = 0.53 from the observed A4481 cores, and 
an excitation temperature of 7200K, we find £f < 0. 

The strong lines in Fig. [5] fall between theoretical curves 
for the parameter a = 7a/(2 • AAd) of 0.01 and 0.1. Note 
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that 7 (without a subscript) traditionally means 7^ , or the 
FWHM of a dispersion profile expressed as a function of 
cu = 2nc/X. 

Using T = 8300K, and a mean A = 4.6 • 10~ 5 cm, we 
obtain 7 = 4.29 • 10 9 s" 1 for a = 0.1. This is 41 times 
the classical damping constant 7^ = 0.22/A 2 s^ 1 (with A 
in cm). For a = 0.01, the result is 4 times the classical. A 
value commonly cited for the empirical damping constant is 
10 times the classical value (cf. Mihalas 1970). Cowley and 
Cowley (1964) obtained 7/701 = 25 for the sun. We thus 
see the value from this crude curve of growth is entirely 
reasonable for a star on or near, the main sequence. 

7. Abundances 

All abundances were determined from depth-dependent 
model calculations of equivalent widths assuming zero mi- 
croturbulence. Oscillator strength references not explicitly 
mentioned are from NIST (Ralchenko, 2010, preferred) or 
VALD (Kupka et al. 1999), except for Zr II, where the val- 
ues are from Ljung, et al. (2006). 

Plots of abundance vs. log(Wx) were made for all species 
with more than three lines. Outliers could nearly always 
be understood and either corrected or dropped because of 
errors, such as misidentifications or blends. Whenever plots 
of abundance vs. equivalent widths exhibited a downward 
trend for stronger lines, abundances were determined from 
lines with equivalent widths < 20 mA. 

Additional plots were routinely made of abundance vs. 
excitation potential, and equivalent width vs. excitation po- 
tential. Our restriction to weak lines only for abundances 
eliminated most systematic effects, though some small ef- 
fects are inevitable. Additional plots, and detailed line-by- 
line results are available from the first author or at the url: 
www. astro. umich.edu/^cowley/hdl01412/ 

Abundances are summarized in Tabled Errors in Col 3 
are standard deviations for the number of lines used, or the 
difference in abundances when only 2 lines were available. 
For V 1, the error is the difference from a single V 1 line, 
and the mean for V II. 

No exotic elements, were identified. Barium is the heav- 
iest element positively identified, and the abundances listed 
for the lanthanides cerium and europium are upper limits. 
There was no indication of gallium or the heavier noble 
gases. 

Logarithmic abundance differences-star minus sun-are 
plotted in Fig. [71 The abundances may reflect a mild A Boo, 
or Vega-like abundance mechanism, where the refractory 
elements are depleted while the most volatile elements are 
more nearly normal (Takeda 2008, see also Adelman, et al. 
2010, in preparation). This is not unexpected. Gray and 
Corbally (1998) classified the Herbig Ae star HD 37411 as 
a A Boo star, while Acke and Waelkens (2004) describe the 
Herbig Ae star HD 100546 as a "clear A Bootis star," based 
on their abundance analysis. 

Neither Vega nor HD 101412 would qualify as a A Boo 
star, as they do not meet the 1 dex iron-peak deficiency of 
Paunzen (2004) for the A Boo stars. Our average deficiency 
for the elements with 50%Tc > 900K is about 0.5 dex, 
half the depletion for classical A Boo stars. Nor is there an 
obvious trend within these refractory elements with 50%Tc. 

Our suggestion that the abundance pattern resembles 
that of the A Boo stars rests on average abundances for the 



Table 2. Adopted abundances. Upper limits assume a 
1 mA line for Ce 11 A4686, and Eu 11 A4205. The 50% 
Condensation temperatures (50%Xc) in the table are from 
Lodders (2003). Elements with low values of 50%Tc are 
volatile. Solar abundances are from Asplund, et al. (2009). 



Spec 


lout El /Ntnt) 

iw O\"^ / J " TOT j 


sd. 


No. 


Sun 


50%Tb(K) 


C 1 


—3.70 


0.26 


10 


—3.61 


40 


N 1 


—3.49 


0.26 


10 


—4.21 


123 


O 1 


—3.21 


0.19 


11 


—3.35 


179 


Na 1 


—5.94 


0.09 


4 


—5.80 


953 


Mff T 
1VA 6 1 


—5.00 


0.24 


11 


—4.44 


1327 


Mff TI 
1VA & 11 


—5.06 


0.19 


8 






Al 1 


—6.22 


0.56 


5 


—5.59 


1641 


Si 1 


—4.92 


0.40 


9 


—4.53 


1302 


Si 11 


—5.30 


0.31 


10 






S 1 


—4.96 


0.10 


8 


—4.92 


655 


Ca 1 


—6.38 


0.07 


3 


—5.70 


1505 


Ca 11 


—6.30 


0.52 


10 






Sc 11 


—9.36 


0.09 


11 


—8.98 


1647 


Ti 11 


—7.62 


0.16 


43 


—7.09 


1573 


V 1 


—8.29 


0.19 


1 


—8.11 


1427 


V 11 


—8.48 


0.22 


28 






Cr I 


—6.89 


0.26 


17 


—6.40 


1427 


Cr 11 


—6.95 


0.16 


6 






Mn 1 


—6.98 


0.32 


12 


—6.61 


1150 


Mn 11 


—6.95 


0.13 


7 






Fe 1 


—5.07 


0.19 


18 


—4.44 


1328 


Fe 11 


-5.03 


0.17 


30 






Co 1 


-7.62 


0.17 


17 


-7.05 


1347 


Co 11 


-7.46 


0.24 


6 






Ni 1 


-6.34 


0.21 


40 


-5.82 


1348 


Ni 11 


-6.12 


0.27 


2 






Zn 1 


-8.54 


0.12 


3 


-7.48 


723 


Sr 11 


-9.27 


0.36 


2 


-9.17 


1455 


Y 11 


-10.34 


0.19 


9 


-9.38 


1647 


Zr 11 


-9.73 


0.24 


16 


-9.46 


1736 


Ba 11 


-10.68 


0.09 


2 


-9.86 


1447 


Ce 11 


< -11.22 




1 


-10.46 


1477 


Eu 11 


< -12.37 




1 


-11.52 


1347 




II I I I I I I I I I I I I I I I I I u 

500 1000 1500 

50% Condensation T(K) 

Fig. 7. Logarithmic abundance differences vs. 50%X"c- 
Error bars are from Table [2] Volatile elements, C, N, O 
are enclosed within the ellipse on the left; refractory ele- 
ments fall within the ellipse on the right. The intermediate 
volatiles, Zn and S, are between the dashed vertical lines. 
The decided outlier, Zn, is discussed in the text. 
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elements with volatility extremes. The centroids of the el- 
lipses in Fig.[7]are significantly displaced from one another, 
indicating the most volatile elements are not depleted while 
the most refractory elements are. 

7.1. Outliers 

The abundance pattern of Fig. [7] cannot be explained in 
terms of condensation temperature alone. 

7.1.1. Nitrogen 

The high nitrogen abundance may be partially due to 
NLTE. Kamp, et al. (2001) find nitrogen quite variable 
among A Boo stars, and subject to NLTE. For example, 
the LTE abundance excess of N for HD 75654 was +0.65 
dex. The NLTE calculation reduced this value to +0.30 dex. 

7.1.2. Sodium 

Sodium (50%Tc = 953K), might be considered an interme- 
diate volatile. We find its depletion in HD 101412 to be only 
0.14 dex. However, Paunzen, et al. (2002) reported sodium 
variations from —1.3 to +1.2 dex with respect to the sun. 
Kamp et al. (2008) discuss this further, and note mecha- 
nisms beyond equilibrium condensation that might account 
for the scatter. 

7.1.3. Zinc 

Paunzen(2004) listed zinc among the elements expected to 
be deficient in the A Boo stars. To the extent that we sug- 
gest the abundances in HD 101412 resemble those of the A 
Boo stars, our finding of an depletion of 1.06 dex ceases to 
indicate an outlier, and supports the assertion. At the same 
time, it complicates an interpretation based on volatility. 
Surely factors other than equilibrium condensation affect 
the stellar abundances of zinc as well as sodium. 

8. Conclusions 

The atmospheric structure of HD 101412 deviates from 
a classical model (e.g. ATLAS9). But the deeper photo- 
sphere, as probed by weaker lines, has a line spectrum 
well produced by the theoretical model with T e g 8300K, 
log(g) — 3.8. Elements with lines from the first and second 
spectra yield abundances in good or satisfactory agreement. 
In the worst case, silicon, Si I and Si n disagree by 0.38 dex, 
only a factor of 2.4. 

The atomic lines are also in general agreement with a 
mean excitation temperature in the range 7200 to 7750K. 
An excitation temperature in this range is expected from a 
model with the adopted parameters. 

Strong lines show an anomalous saturation, which we 
have attributed to maximum line depths, ro, less than 
that predicted by the ATLAS9 model. Similar informa- 
tion comes from curves of growth, where even the observed 
depths of the Mg n A4481 doublet (0.53) would lead to a 
imaginary value for the microturbulence. 

Atomic emission features show departures of the upper 
atmosphere from classical. In addition to emissions seen in 
the low Balmer members, broad, generally weak emissions 
appear, for example, in [O i] A6300, Na I Di and D2, and 



the O 1, triplet AA7772, 7774, and 7775. The widths of these 
features are readily measured. The wavelength spread from 
the violet wing to the red wing, of atomic emissions corre- 
spond to velocities of 100 to 200 km s^ 1 . We suggest the 
upper atmosphere is heated by infalling material from a pri- 
mordial disk, and that this heating is probably responsible 
for the anomalous saturation. 

Most abundances are less than 1 dex from solar values. 
With the notable exception of zinc, there is a suggestion 
that refractory elements are depleted. The volatiles are nor- 
mal, or in the case of nitrogen, enhanced. 
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